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ABSTRACT 

A simple integral relation is obtained for the distribution functions of jet speeds and 
jet-counter jet flux asymmetries of classical double radio sources on the basis of a sim- 
ple model in which it is assumed that the jets are intrinsically symmetric and that 
the asymmetries are attributed to relativistic beaming. Analytic expressions relating 
the mean jet speed and the mean jet flux asymmetry, and their variances are derived. 
The results are considered in the light of orientation-based unified schemes, and an 
analytical equation for the critical angle separating quasars and radio galaxies is de- 
rived. The mean jet speeds of classical double radio sources are estimated using the 
asymmetry of jet-counterjet flux densities taken from the 3CRR and B2 samples. For 
FRI radio sources the mean jet speed is ~ (0.54 ± 0.03)c, while for FRII low-redshift 
radio galaxies and intermediate redshift quasars the values found are ~ (0.4 ± 0.06)c 
and <; 0.6c respectively. 
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1 INTRODUCTION 

There is observational evidence that the jets of double FRI 
and II radio sources are initially relativistic. On parsec scales 
this idea is supported by the observation of superluminal 
and relativistic motions, rapid variability and high bright- 
ness temperatures. On kiloparsec-scales the correspondence 
of jet-sidedness with depolarization asymmetry of both FRI 
(Morganti et al. 1997) and FRII radio sources (Laing 1988, 
Garrington et al. 1988) is interpreted as the brighter jet ap- 
proaching the observer, thus supporting the idea of relativis- 
tic beaming of the radiation of the jet. In all double radio 
sources in which one-sided parsec and kiloparesc-scale jets 
are observed, the parsec-scale jet points in the direction of a 
kiloparsec-scale jet, implying that relativistic speeds persist 
on kiloparsec scales (Pearson & Readhead 1988, Wardle and 
Aarons 1997). 

Other observational data are consistent with this pic- 
ture. Laing et al. (1999) analysed the asymmetries in the 
jet flux densities of 38 FRI radio sources selected from 
the B2 sample of Parma et al. (1987) and concluded that 
jets are intrinsically symmetric and initially relativistic with 
(3 ~ 0.9c. Analysing a large sample of FRII radio galax- 
ies with z < 0.3 Hardcastle et al. (1999, hereafter H99) 
found that relativistic beaming is needed to explain the re- 
lationship between core and jet prominences, with speeds 
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between (0.5 — 0.6)c on kiloparsec scales. Wardle & Aaron 
(1997, hereafter W97) analysed the jet-counterjet flux ra- 
tios of 13 3CR FRII quasars imaged by Bridle et al. (1994). 
They showed that only modest intrinsic jet asymmetries are 
allowed and that the large jet-counterjet brightness ratios 
can be attributed almost entirely to the Doppler effect. The 
idea that jets are intrinsically symmetrical is supported by 
other observational evidence, for example, the detection of 
two parsec-scale jets in the giant double radio source DA240 
with a flux density ratio of unity (Saripalli et al. 1997) and 
the existence of the Compact Symmetric Objects which are 
thought to be formed by symmetric jets. Readhead et al. 
(1996) and Owsianik & Conway (1998) argued that Compact 
Symmetric Objects represent a very early phase of classical 
double radio sources. 

These analyses have assumed that the jet speeds associ- 
ated with relativistic beaming can be adequately described 
by a single jet speed for a given source. This seems to be a 
reasonable approximation in view of the remarkable symme- 
try of the hot spots observed in the extended lobes of FRII 
sources, which need to be constantly supplied with energy 
by such jets over the lifetime of the source. There is how- 
ever an issue whether or not the observed emission from the 
most compact jets is associated with a single jet speed. The 
polarisation structures of the compact jets observed in FRI 
sources by Laing (1993), Hardcastle et al (1996) and Laing 
and Bridle (2002) suggest that these jets can be modelled by 
a spine-shear layer structure, in which entrained material in 
the shear layer would result in a dispersion in the speeds of 
the radio emitting material. These structures are assumed to 
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be fast decelerating jets which do not persist much beyond 
their observed locations in the radio galaxies. These effects 
become important beyond what they term the flaring point 
of the jets. Most of the emission from the compact jet origi- 
nates, however, from the jet spine which is assumed to have 
a constant speed. Thus, even in the case of the FRI jets, it 
would be an adequate assumption that the bulk of the emis- 
sion is associated with the jet spine, certainly at distances 
less than the flaring point. It will be assumed that jet speeds 
within a given FRII source can be adequately modelled by 
a single jet speed in the analysis which follows. 

To date, estimates of the jet speeds have been found 
from Monte-Carlo simulations in which the jet properties are 
determined by relativistic beaming and consistency sought 
between the data and the model. In this paper, the recon- 
struction of the distribution function of jet speeds is found 
using an inverse-problem approach, in which the mean jet 
speed and its standard deviation are found directly from 
the observed jet-counterjet brightness distribution of FRI 
(Laing et al. 1999, hereafter L99) and FRII radio sources 
(Bridle et al. 1994, Hardcastle et al. 1998, hereafter B94 
and H98). In Section 2, the inverse problem is formulated 
and the distribution function of jet speeds and analytical 
equations for estimating the mean jet speed and its variance 
are derived. These results are examined in the context of 
orientation-based unified schemes in Section 3, and a new 
method of estimating the critical angle which separates ra- 
dio galaxies and radio quasars is derived. In Section 4 and 5 
the reduced equations are applied to the data available for 
FRI and FRII radio sources, and the discussion and con- 
clusions are presented in Section 6. In our computations, 
Hubble's constant is taken to be Ho = 50 km s _1 Mpc -1 
and the world model Qq = 1, =0. 



2 THE INVERSE PROBLEM APPROACH 
2.1 The model 

Although superluminal motions have been interpreted as 
"pattern speeds" rather than the bulk motions of the jets 
themselves, we adopt the latter picture in the present analy- 
sis as the most natural interpretation of the observed asym- 
metries of the jets since it is entirely consistent with other 
constraints on viable models for the physics of FRII sources. 
The specific assumptions made in the model are as follows. 
First, it is assumed that, intrinsically, the jets are symmet- 
ric and that the asymmetric observed brightness distribu- 
tion is entirely attributable to relativistic beaming. Second, 
it is assumed that jets are oriented randomly with respect 
to the line of sight. The primary selection criterion for the 
radio sources is that they are chosen from complete flux- 
density limited samples of extended radio sources at a low 
radio frequency of 178 MHz. Since the flux densities of the 
sources are wholly determined at this frequency by the to- 
tal flux density of the unbeamed extended emission with 
no contribution from Doppler-boosted emission, there is no 
bias in the selection of the radio jets. Third, we can adopt 
the procedures we have already discussed to understand the 
asymmetries of distributions of hot-spots to analyse the sta- 
tistical properties of the flux density ratios of the radio jets 
(Arshakian and Longair 2000). 



2.2 The distribution function of jet speeds 

Assuming that the jets are intrinsically symmetric and that 
the observed jet flux asymmetries are entirely due to Doppler 
beaming, the ratio of the jet and counterjet flux densities, 
J = Sj/Scj, can be written 

where is /3jC is the speed of the jet, which is inclined at 
an angle 9 £ [0, 7r/2] to the line of sight to the observer, 
and 5 = m + a; the constant m = 2 for a continuous jet 
(Scheuer & Readhead 1979). a is the spectral index, defined 
by S v oc v~ a and taken to be 0.6. 

Let us invert equation (1) to express to — f3jcc*s9 in 
terms of J, 



jl/6 



/3j cos 9, 



(2) 



" JV« + 1 

Suppose the probability distribution of jet flux asymmetries, 
as parameterised by to, is given by the distribution function 
g(w) and is known from observation. Our task is to deter- 
mine the distribution function of jet speeds G(/3j) given g(u>). 

According to the model, /3j and cosf? are independent 
random variables since there is no bias in the selection of 
the sources and we assume that the jet can be characterised 
by a single jet speed. Therefore, the cumulative distribution 
function D(ui) = g(x) dx for lo is, 



D(u>)= J J G(P i )K(9)d[3 i d9. 



(3) 



Assuming the radio axes are distributed isotropically over 
the sky, K{9) = sin 9. Opening the limits of integration, 



■Jui Pi 

and differentiating with respect to w, 
G(w) = -wg'(u). 



(4) 



(5) 



This procedure for relating G{u>) and the derivative of g(u>) 
is similar to the reduction carried out by Arshakian and 
Longair (2000). 

2.3 The mean jet speed and variance 

Multiplying equation (4) by oj k du and integrating from zero 
to 1, we obtain the k th moment of the distribution g(u>), 



Vk 



fc + 1' 



(6) 



where uok is the fc th moment of the distribution function of 
true velocities, 



H>* = / 1 ^G(/3 j )d/3 j . 
Jo 



(7) 



The first and second moments of the distribution function 
for the jet speeds are equal to the mean speed and the mean 
square speed of the jet, i>oi = «o and V02 = Pf, while the 
same moments for the function g(u>) are vi = uJ and V2 = u> 2 . 
Hence, from equation (6), the mean jet speed is given by 



Pi = 2aJ, 



(8) 
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while the mean squared jet speed is (3f = 3ui 2 . The disper- 
sion of jet speeds is given by, 



4 3 =3^_4U 2 . 



(9) 



3 UNIFIED SCHEMES 

Scheuer (1987) suggested that powerful radio galaxies and 
radio quasars arc drawn from the same distribution of 
sources, but are viewed at different angles, the quasars hav- 
ing smaller angles to the line of sight than the radio galax- 
ies. There is considerable evidence supporting this idea from 
analyses of the average sizes, structural asymmetries, mis- 
alignment angles and the jet-sidedness of radio quasars and 
radio galaxies (Barthel 1989, Scheuer 1995, Best et al. 1995, 
Arshakian & Longair 2000). Let us derive the distribution 
function of jet speeds and the relations between the mo- 
ments of the distribution functions of jet speeds and g(uj) 
within the context of this orientation-based unified scheme. 



3.1 General case 

Suppose the radio axes of certain FRII sources are observed 
in the range of angles 9 6 (8i, O2) and that the distribution of 
angles is random in this interval. Then, K(6) d9 = C sin 6 d6, 
where C = l/(cos#i — cos 62). Then equation (4) becomes 



g(u) = C { 



o j cos $2 



o j cos 6 
1 

/ 

k uj j cos Q\ 



Pi 

G(ft) 

Pi 



d/3j , for uj < cos 62 ; 



d/3j , for it) > cos 62 , 



(10) 



under the corresponding conditions that /3j £ [0, 1], cosf? £ 
[cos 61, cos 62] and uj 6 [0, cos#i]. 

This system of equations was obtained independently 
by Baryshev and Teerikorpi (1995) in their study of the 
distribution function of hotspot velocities for FRII radio 



where C'g = 1/ cos# c . For quasars, an analytical solution 
for the moments can be obtained from the first integral in 
equations (11), 



_ (n+lXl-cosflc) 



(13) 



Hence the mean speed, mean squared speed and the disper- 
sion are, 



_ - 2(l-cosfl c ) 

V01 = (Pj)Q = 9a HQ, 



1 — cos 2 



.--2 _ 3 (1 - cosflc) — -. 



„02 EE (#)q = ^_ 



cos 3 6 C 



(^ 2 )Q, 



and, 



K)Q = (/3f)Q-(/^)Q. 



(14) 
(15) 

(16) 



respectively. From equation (12), the corresponding mo- 
ments for radio galaxies are 



VQ n 



n + 1 

cos n e c 



and so, 



(A)g = 



COS t 

3 

cos 2 8, 



■ (^)g, 



and, 



(17) 

(18) 
(19) 

(20) 



If radio galaxies and radio quasars are intrinsically 
the same objects, the mean jet speed of radio galaxies 
and quasars drawn from the same sample should be equal: 
(/3j)q + g = (/3j)g- From equations (8) and (18), 



(50c 



.(w)q+g. 



(21) 



3.2 The mean jet speed and variance 



3.3 Test of analytical equations 



In the orientation-based unification scheme, the radio galax- 
ies and radio quasars are observed at different angles to the 
line of sight and are separated by a critical angle 6 C ~ 45° 
(Barthel 1989). For quasars 6\ = and 62 = C , and so 
equations (10) become, 



J j cos i 



?H = C Q < 



1 

/ 



Pi 

G(/3j) 
Pi 



d/3j , for lo < cos 9 C 



(11) 



d/3j , for ui > cos ( 



where Cq = 1/(1 — cosd c ). For the radio galaxies 6 1 = 6 C , 
62 — tt/2 and so, 



g(u>) = C G 



G(Pj] 
Pi 



d/3j , for ui < cos 9 C 



(12) 



uj j cos 6 C 
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In order to test the formalism involving equations (8,9), 
(14,16), (18,20) and (21), we have modelled the distribu- 
tion function g(ui) of jet flux asymmetries from equation (2) 
by generating 5,000 random values of /3j and under the as- 
sumptions that (i) the distribution function of jet speeds is 
Gaussian with mean /3 m and standard deviation ap m and (ii) 
cos# is uniformly distributed in the ranges [0, 1] for equa- 
tions (8,9), in the range [0,cos# c ] for equations (14,16) and 
in the range [cos 9 C , n/2] for equations (18,20). For various 
combinations of the free parameters of a model, /3 m , cr/3 m and 
# c , the distribution function of uj has been modelled for (i) a 
joint sample of radio galaxies and quasars, (ii) radio galaxies 
and (iii) quasars. The mean and mean square values of g(u>) 
for these samples are used to calculate the mean jet speed, 
standard deviation and critical angle using the above ana- 
lytic equations. We find very good agreement between the 
different combinations of free parameters and the calculated 
mean jet speed, standard deviation and critical angle. 
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4 JET SPEEDS IN FRII SOURCES 

In order to apply the results of Sect. 3, high-resolution, high- 
sensitivity maps of complete samples of sources are needed. 
The samples of sources used to study the FRII radio sources 
are listed in Table 1, which provides statistical data on the 
total numbers of sources and the numbers of different types 
of source within the samples. The two samples consist of a 
sample of high-redshift 3CR FRII quasars studied by B94 
and a sample of low-redshift radio galaxies with z < 0.3 
studied by H98. These studies adopted the same definition of 
what is meant by a jet and the same method of determining 
the jet flux densities which are used in the following analysis. 

The problem in analysing these data is that the detec- 
tion of jets is limited by the sensitivity and resolution of the 
observations. In only a few cases are measured values of J 
available. The following criteria have been adopted. We con- 
sider only sources with definite (D), definite/definite (DD), 
definite/probable (DP) and probable/probable (PP) jets de- 
fined for the whole jet length, following the terminology of 
H98. The sources for which there are only upper limits to 
the straight jet and counterjet fluxes are classified as 'inde- 
terminate' and those in which the straight jet flux density 
of the brighter jet is less than the straight counterjet flux 
density were 'excluded'. We take the brighter jet to be that 
approaching the observer and use flux densities from the 
straight part of the jet for estimating the J and hence w. 
The statistics of the jets in the sample are given in Table 1. 



4.1 FRII radio galaxies 

We analyse the jet velocities of the 32 high-excitation and 12 
low-excitation radio galaxies separately, because it has been 
shown that the FRII low-excitation radio galaxies form a 
separate group of radio galaxies which are not unified with 
high-excitation sources (Laing et al. 1994, Arshakian & Lon- 
gair 2000). 

H98 state that their low-redshift sample of 44 radio 
galaxies is free from orientation bias, as the 3CR sample 
was selected on the basis of a low-frequency (178MHz) sur- 
vey in which the jets contribute a negligible fraction of the 
total flux density. The low redshift limit [z < 0.3) means 
that there are no quasars in the complete sample and so the 
effects of orientation are associated with a sample of radio 
galaxies alone. Of the 32 high-excitation sources, estimates 
of J could be made for 17 sources - fluxes were measured 
from both jets for 5 sources, and only these have accurate 
values of lo, while for the remaining 12 only lower limits are 
available. Throughout the paper we treat the lower limits 
as representing the true values and so this will lead to an 
underestimate of the mean velocity. Most of the 15 indeter- 
minate sources are assumed to be inclined at a relatively 
small angle to the plane of the sky, and are not observable 
because their flux densities are not Doppler-boosted. There 
is therefore a lack of sources at large angles to the line of 
sight among the 17 for which J values can be estimated. If 
we ignore this bias and assume that the 17 radio galaxies 
have an isotropic distribution, then an upper limit to the 
mean jet speed can be estimated from equation (8), since 
values of J close to one are under-represented. The mean 
value of u> for the 17 sources is log = 0.24 and the mean jet 
speed is < 0.49c. 



Let us estimate the lower limit to the mean jet speed 
by considering only those 16 sources of the 32 with one def- 
inite jet or definite/probable jets defined by the whole jet 
length. We have argued that these 16 sources are inclined 
preferentially towards the line of sight to the observer and so 
we assume that they are viewed within a maximum viewing 
angle 6 C to the line of sight, while the remaining 16 radio 
galaxies have orientations between 9 C and ir/2. An estimate 
of the critical angle then allows the mean jet speed and the 
standard deviation to be calculated from equations (14, 16). 
On the assumption that the sources with definite and defi- 
nite/probable jets are oriented isotropically between 0° and 
9 C , we find, 



1 



N B + Nop 



60°, 



(22) 



where No and Nop are the number of sources with definite 
or definite/probable jets respectively. 

For the 16 high-excitation radio galaxies with D/DP 
jets, the value of J can be estimated only for 11 sources for 
which the mean value of lo is log = 0.29, and so the mean jet 
speed on kiloparsec scales is f3 G ~ 0.4 ±0.06 from equations 
(14) and (16) if 6 C = 60°. The mean jet speed lies in the 
range (0.35 - 0.43)c ± (0.06)c if 6 C G [50°, 70°]. From this 
analysis it may be concluded that the mean jet speed of 
low-redshift FRII radio galaxies is about 0.4c with an upper 
limit of 0.5c. This range of speeds is slightly lower than those 
obtained from the analysis of jet prominences by H99, who 
found v m (0.5 — 0.7)c. Despite the number of upper limits 
in their analysis, they believe that these values are not far 
from the true ones. 

A similar analysis can be carried out for the sample of 
8 definite and 5 possible jets in low-excitation radio galaxies 
to estimate mean and upper limits to the mean jet speed. 
The source 3C15 has an excessively large value of lo = 0.67, 
which is not typical of two-sided jets. Excluding it from the 
final sample and using equations (8) and (14), the 7 sources 
result in a probable mean jet speed, v <; 0.26c, and an upper 
limit, v 5s 0.41c, if 6 C ~ 55° for these D/DP jets, assum- 
ing that low-excitation radio galaxies form an independent 
group of sources with a uniform distribution of orientation 
of the radio axes on the sky. This result suggests that on the 
average the jet speeds of low-excitation galaxies are slightly 
smaller than those of high-excitation galaxies, but since it 
is based on small number statistics, it is only a provisional 
conclusion. 

Arshakian & Longair (2000) have shown that the lobes 
of low-excitation radio galaxies have non-relativistic expan- 
sion speeds, which are less than the speeds of narrow-/broad- 
line radio galaxies at the same redshifts. This result is consis- 
tent with the lower jet speeds of low-excitation radio galaxies 
and/or the tendency for them to belong to denser environ- 
ments (Hardcastle & Worral 1999) which may be the reason 
for lower expansion speeds of radio lobes. 



4.2 FRII radio quasars 

The case of the quasars is more complicated. Quasars with 
extended radio structures were selected by B94 for VLA ob- 
servations at 5 GHz and so their radio axes must lie at some 
angle greater than # m i n to the line of sight. This inevitably 
leads to the selection of quasars oriented at angles to the line 
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Table 1. The properties of the radio sources in the B94 and H98 samples 



Properties 


3CR quasars (B94) 


Low Rcdshift Radio Galaxies with z < 0.3 (H98) 


Total number of sources in sample 


13 


44 


High excitation radio galaxies Low excitation radio galaxies 


Number of sources in sub-sample 


13 


32 12 


Number of sources with detectable jets 


13 


22 8 


Sources for which J can be estimated 


12 


17 7 


Sources for which J is indeterminate 




15 


Excluded sources 


1 


5 1 


Sources for which J can be estimated 




Broad-linc Narrow-line 



5 12 



of sight closer to the upper limit at which they would still 
be classified as radio quasars, as compared with an unbiased 
sample - the sources lie in the range 6 £ [# m in, C ]. For a crit- 
ical angle in the range 45° to 60°, W97 estimate the most 
probable range of jet speeds to be (0.6 — 0.7)c for 13 quasars 
observed by B94 by considering the lower limits to the jet- 
countcrjet brightness ratios to represent their true values. 
For the same sample, we calculate from equations (14) and 
(16) the range of mean jet speeds to be ~ (0.5 — 0.56)c and 
<Tf3 Q ~ 0.22, if 45° < 6 C < 60°. This result supports the esti- 
mates of W97, confirming the presence of highly relativistic 
speeds in the jets of FRII quasars. 

Among the sources in this sample of quasars, analysis 
of the jet-counterjet flux asymmetries and the luminosity of 
FRII quasars shows that, for at least two quasars (3C68.1 
and 3C208), the jet flux asymmetry cannot be entirely at- 
tributed to Doppler beaming and contains a significant frac- 
tion of confusing emission from the straight part of the jets 
(see Appendix). This leads to smaller values of u>, which 
results in an underestimate of jet speeds. Excluding these 
quasars from the final sample results in a total of 11 FRII 
quasars. The lower and upper limits to the mean jet speeds 
can be estimated assuming that the quasars are inclined 
preferentially, (i) at small angles (0 — 6 C ) and (ii) at large 
(#min — 90°) angles, using equations (14) and (18) respec- 
tively. Taking ZJ ~ 0.47 for all 11 FRII quasars, we find the 
upper limit to be /3q > 1 for # m in > 20° from equation (18), 
showing directly that the radio quasars cannot be oriented 
at large angles to the line of sight - if the quasars were ori- 
ented at large angles to the line of sight, there should be 
many more quasars with double-sided jets with small values 

Of LO. 

Several analyses have shown that 8 C — 45° for FRII ra- 
dio sources in the redshift range from 0.5 — 1. 8 of the 11 
quasars in our sample lie in the redshift range 0.4 — 1.1, and 
so we can adopt the critical angle (45°) as an upper limit to 
their inclination angles. From equations (14) and (16), the 
mean value of lo ~ 0.47 results in a mean jet speed of be- 
tween (0.55c - 0.63) ± 0.05c, if 9 C = (45° - 60°), which is in 
fact a lower limit to the mean jet speed. This result demon- 
strates that the jet speeds on kiloparsec-scales increase from 
low-redshift FRII radio galaxies, v <i 0.5c, to intermediate- 
redshift FRII quasars, v <; 0.6c. 

The speed v ~ 0.8c is obtained by considering the 
jet flux asymmetries of 5 low-luminosity quasars, 3C 215, 



249.1, 334, 336, 351 which all possess two-sided jets and 
P < 10 28,8 WHz~ 1 . These quasars are probably oriented 
close to a critical angle, 45°, or else two sided jets would 
scarcely be detectable. Then, (3q = 0.56/ cos 45° w 0.8. 

5 JET SPEEDS IN FRI RADIO GALAXIES 

L99 conducted a detailed investigation of 38 weak FRI ra- 
dio galaxies selected from the B2 sample of Parma et al. 
(1987). They defined a sample free from orientation bias by 
restricting attention to those sources with total flux densi- 
ties, excluding the radio core, exceeding the survey limit at 
408 MHz. All these sources were known to have well-defined 
straight jets satisfying the criteria of Bridle & Perley (1984). 

Assuming that jets are intrinsically symmetric, they 
analysed the jet-counterjet brightness ratios at the flaring 
point of the main jet to estimate the range of jet veloci- 
ties. The lower limits of the brightness ratio for 10/38 FRI 
sources were treated as measured values. They analysed the 
distribution of the logarithm of jet-counterjet flux ratios by 
fitting it to models using a maximum likelihood approach. 
Single-velocity models were found to be unsatisfactory, but 
a best-fitting velocity was found to be /3r = 0.72. For their 
spine/shear- layer model, the common features are a maxi- 
mum velocity ~ 0.9c and a large dispersion <; 0.6c. They 
found that the range of velocities at the flaring point is 
0.56 </3a< 0.8. 

For the same sample, we have calculated the jet prop- 
erties using equations (8) and (9) assuming that the jets are 
uniform at the flaring point, that is, our analysis is indepen- 
dent of the jet model. We use the jet flux asymmetry pa- 
rameter u) defined by (2) with 8 = 2.6. The mean jet speed 
at the flaring point was found to be cf3 a — (0.54 ± 0.03)c 
and the standard deviation of jet speeds cr^g = 0.18c. The 
higher moments of the distribution determine the skewness 
and kurtosis of the distribution function of jet speeds, which 
are Skewness = /?|/(/?£) 3/2 = 4^/(3^2) 3/2 = 0.64 and 
Kurtosis = ]4/(f4) 2 - 3 = 5Z7/(3^) 2 - 3 = -1.7. The 
negative kurtosis indicates that the distribution function of 
jet speeds is flat-topped, while the positive skewness im- 
plies that the distribution is asymmetric with the maximum 
shifted to slightly smaller velocities. The maximum speed, 
v <; 0.72c, is found for the source 1652+39. Despite the fact 
that the mean jet speed is within 3cr of the high range of 
speeds 0.56c — 0.8c obtained by L99, we conclude that most 
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FRI radio sources have jet speeds less than 0.54c, with a 
large spread from subrelativistic to high relativistic speeds. 

L99 also point out a selection effect leading to a lack of 
face-on sources because their jets may be too faint to detect 
or to measure accurately. Inclusion of these sources would 
decrease the mean jet speed even more, [3 a < 0.54. 



6 CONCLUSIONS AND DISCUSSIONS 

Subject to the simplifying assumptions of the model, which 
are described in the first paragraph of Section 2.1, the prin- 
cipal results of this analysis are: 

(i) On the basis of an intrinsically-symmetric relativistic 
jet model, the distribution function of jet speeds is derived, 
and analytic expressions for the mean jet speed and its stan- 
dard deviation are obtained. The equations for estimating 
the same quantities and the critical angle are also derived 
in the context of orientation-based unification schemes. 

(ii) Applying this analysis to the jet flux asymmetries of 
double radio sources, we find that: 

(a) the mean jet speed of low-redshift FRII radio galaxies 
lies in the range ~ (0.4c — 0.5c) ± 0.05c and the standard 
deviation ~ 0.15c; these speeds are in agreement with 
the range of speeds obtained from the analysis of jet 
prominences, 0.5c — 0.7c (H99); 

(b) a lower limit to the mean jet speed of intermediate- 
redshift FRII quasars is v <; 0.6c, which supports the 
presence of highly relativistic jets in quasars (B94, W97), 
showing that the typical jet speed increases with redshift 
and /or radio luminosity; 

(c) the mean jet speed of FRI radio galaxies at the flaring 
point is about 0.54c with standard deviation 0.18c; the 
peak of the distribution function of jet speeds is shifted 
to lower speeds. 

When the lower limits to cu are taken as exact values, the 
analytic methods of the estimation of the mean jet speeds of 
FRI/FRII sources give values slightly lower than those based 
on testing (W97, H99) or fitting (L99) models to the data. 
In the case of low-redshift radio galaxies, this difference can 
be attributed to the different asymmetry parameters used - 
these are the jet/counterjet flux asymmetry and the jet/lobe 
emission contrast (H99). Despite this, the calculated mean 
jet speeds are in agreement within the limits of errors. Jet 
flux asymmetries have been used in analysing the jet speeds 
in FRII quasars and FRI radio galaxies. Our estimates of 
/3j are consistent with the values of W97 and only slightly 
different from the values of L99. 

An important issue is the following: what are the real 
values of the jet flux asymmetries as parameterised by u>? 
FRI and FRII low-redshift radio galaxies are seen at angles 
between 0° — 90°, while most FRII quasars are inclined at 
angles less than 45° . Relativistic boosting of the jet emission 
is more important for quasars (B94) than for radio galaxies 
(H98), and hence a difference between the lower limits to 
the jet flux asymmetry and the real values should on average 
be higher for quasars than for radio galaxies. Simulating the 
effects of deeper observations, we multiply the lower limits to 
the jet-counterjet flux ratios of FRI and FRII radio galaxies 
by 1.5 and estimate the mean speed in the range ~ (0.6 ± 
0.03)c and ~ (0.5c-0.6c)±0.05c. For FRII quasars, when the 



limits on the flux ratios are doubled, a mean speed between 
0.7c and 0.8c is more acceptable. 
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Figure 1. (a) Jet and (b) counterjet flux densities of 13 quasars versus the luminosity. Filled asterisks are quasars with two sided jets, 
while unfilled asterisks are quasars with one sided jets. Upper limits to the straight jet flux are denoted by J.. 



APPENDIX: EVIDENCE FOR CONFUSING 
EMISSION IN THE COUNTERJET 
BRIGHTNESSES OF QUASARS 

There is a negative correlation between the straight jet 
flux and luminosity, Sj oc p-os4±o.i8 ( Fig la ) ne should 
expect that on the average the counterjet flux densities 
also should correlate negatively with luminosity with some 
spread depending on the jet speed and jet orientation, 



S ei =S i ( ] + % COB0 X t <xP~ 0M . 



v 1 — /3j cos ( 



(23) 



In fact no correlation is observed as can be seen in Fig. lb. A 
plausible explanation is that some of upper limits to coun- 
terjet fluxes densities are overestimated. The expected cor- 
relation is shown in Fig. lb by the dashed line, which has the 
slope -0.54 and passes through the point with coordinates 
log S c j = —0.42 and log P = 28.08, the mean values calcu- 
lated for four low-luminosity quasars (P < 10 28 ' 3 WHz -1 ) 
with two-sided jets (Fig lb). It can be seen that the upper 
limits to the counterjet brightnesses are overestimated for 
two quasars, 3C68.1 and 3C208. A plausible interpretation 
is the presence of confusing emission (B94) in the counterjet 
lobe, which leads to an overestimation of the counterjet flux 
density and which in turn seriously underestimates the jet 
flux asymmetry for these quasars. 
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